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1 Introduction 



The high energy scattering in QCD can be most efficiently addressed in terms of colour dipole 
degrees of freedom. In this approach, originally proposed by Mueller pQ, one considers a fast 
moving projectile as a bunch of dipoles created via a classical branching process. High energy 
evolution of the projectile wavefunction can be then found from the QCD generating functional 
P], which obeys a linear functional evolution equation |2]. Though the evolution of the projectile 
wavefunction is independent of target, the interaction amplitude depends on it. 

The energy evolution of the interaction amplitude will be the central object studied in the 
present paper. It was shown by Kovchegov [3j , that if one assumes that all the dipoles produced via 
evolution interact with the target independently of each other, then the scattering amplitude obeys 
a nonlinear evolution equation presently referred to as the Balitsky-Kovchegov (BK) equation 
IHIHj. The BK equation has been the core of numerous analytical [Sj and numerical [Hj studies. 
It has also been applied for phenomenology [7]. 

It was realised in Ref. [2], that the assumption regarding independent interactions can be 
incorrect for many targets, and in particular for realistic nuclei which are not very dense. Target 
correlations are needed in order to account for more realistic structure of targets. It was shown 
in Ref. j2] that the amplitude obeys a linear functional evolution equation, which is the most 
appropriate tool for the introduction of target correlations. In general, the linear functional 
equation for the amplitude with target correlations cannot be represented as an ordinary equation 
and preserves its functional form. Our goal here is to show that the linear functional evolution 
equation for the interaction amplitude for a generic target with correlations is equivalent to an 
infinite chain of linear hierarchy equations. This chain is the hierarchy of Balitsky 0] reduced 
to dipole operators only ^. Thus we clarify the connection between the Wilson loop approach of 
Balitsky and the Mueller dipole model. We also demonstrate that the large A*";, limit, which is 
in the foundation of the dipole model, is insufficient to guarantee a close form equation for the 
scattering amplitude. This new derivation of the Balitsky's chain should be view as our main 
result. 

We also demonstrate a possibility to rigorously derive an evolution for the amplitude without 
restricting ourselves to a single dipole as a specific choice of projectile. Our results are valid for 
any projectile whose evolution is not affected by high density effects. 

At the end of the paper we discuss an ansatz for target correlations which allows one to obtain 
a certain generalisation of the BK equation. This part is very much along the lines of Ref. [2]. 
Similar result has been recently reported by Janik [9 . 

One of our central goals in this letter is to draw more attention to Ref. [2, in which the linear 
functional formalism was developed as well as many other results appearing in this reference. 

^The true Balitsky's hierarchy, as well as its equivalent, the JIMWLK equation 0, also contain non-dipole 
operators, which account for Nc corrections. In the present paper we have nothing to say about Nc corrections 
which will be systematically ignored. 
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2 Balitsky's chain from the QCD generating functional 



We first consider a generic fast moving projectile wliose wave function can be expanded in a 
dipole basis. Contrary to many previous studies we do not restrict ourselves to a single dipole 
as a projectile. We further assume that non-linear effects associated with high dipole densities in 
the projectile wave function can be ignored. This is a strong assumption which eliminates effects 
which might be associated with pomeron loops, so called enhanced diagrams. 

Let us define a probability density P„ to find n dipoles with coordinates ri, 61, r2, 62, ■ ■ - i^iybi, 
. . . rn, hn and rapidity Y in the projectile wave function, and hi denote the dipole's size and 
impact parameter, both are two dimensional vectors. We define P„ as a dimensionfuU quantity. 
The equation for P„ obeys the classical branching process: 

^0; ri,hi, r2,b2, . ..rn,bn) 

{r-i + r„),6i„,...r„_i,6„) (2.1) 

with bin = bi + rn/2 = bn — Ti/2 being imposed via (5-functions. Two terms of Eq. (|2.H) have a 
very simple meaning: the first one describes the decrease in probability to find n dipoles due to a 
decay of one dipole into two of arbitrary sizes. This probability is equal to 

with p being some infrared cutoff, and 0.3 = asNc/ir. The second term shows the increase in 
probability to find n dipoles due to a creation of a new dipole from n — 1 dipoles with probability 

as (ri + ra)^ 

Eq. 1)2.11) has to be supplemented by initial conditions at Y = Yq, specifying a dipole distribution 
in the projectile. In writing the equation 1)2.11) we explicitly used our assumption that there are 
no nonlinear effects in the projectile. Otherwise, we would need to include dipole recombination 
processes. 

The hierarchy ()2.1|1 can be resolved by introducing a generating functional Z 

, n 

Z(Y-Yo; [u]) = X! / Pn{y - Yq] ri,bi, r2,b2,. . . ,ri,bi,. . . ,rn,bn) Yl u{ri,bi) d'^ Vid'^ k 

n=l i=l 

(2.2) 

where u{ri,bi) = Ui is an arbitrary function of rj and bi. It follows immediately from 1)2.1)) that 
the functional ()2.2)) obeys the condition: at u = 1 

Z{Y - Fo; [u = l]) = 1. (2.3) 



dPniY - Yq] ri,&i, r2,&2, ■ ■■rn.br. 
a.dY 



- ^^(r,)P„(r - 
1=1 



+ E fei^ Pn-i - ri, 61, r2, 62, 
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The physical meaning of (|2.3p is that the sum over all probabilities is one. 

Multiplying Eq. ()2.1|1 by the product HiLi integrating over all rj and bi, we obtain the 

following linear equation for the generating functional: 



dZ 

asdY 

with the definitions 
and 



J £rd% V^i^i(r, 6, [u]) Z + j £ r £ r' d% Vi^2{r, r\ b, [u]) Z. (2.4) 



Vi^i{r, b, [u]) = as uj{r) u{r, b) j^jj^ (2-5) 



Vi^2ir, r', b, [u]) = — u{r', b + — - — ) u{r - r',b - -) — — - . (2.6) 

2iT r'^[r — r)^ 2 2 ou{r, b) 

The functional derivative with respect to u{r,b), plays the role of an annihilation operator for a 
dipole of the size r, at the impact parameter b. The multiplication by u{r, b) corresponds to a 
creation operator for this dipole. Eq. ()2.4|1 has an extra b integration compared to the one first 
derived in Ref. |2j. Motivated by the fact that evolution kernels do not depend on the impact 
parameter, we ignored the 6-dependence in Ref. j2j| in order to simplify the presentation. This 
procedure also corresponds to an approximation in which the impact parameter is considered to 
be much larger than any dipole size. The correct analysis has to preserve the true kinematics and 
this amounts to tracing the 6-dependence which is done in Eq. fl2.4|) . Eq. ()2.4|) can be also obtained 
starting from the JIMWLK equation [9j. 

The n-dipole densities in the projectile p^(ri, 6i, . . . , r„, 6„) are defined as follows: 

p^(ri,6i ... ,r„,6„;r - Y,) = - T{ ^ Z {¥ - Y^- [u])\u=i 

n\ f-Ji dui 

One can recast the chain of equations ()2.1|) into hierarchy for p^: 

dff„{ri,bi ... ,rn,K)_ ^ | _ ^ tu(r,) p^(ri, 6i . . . , r„, 6„) + 

\ 1=1 



asdY 

~^ J 271 rf (r; — rr ~i (2 7i]rfri 



i=l 



(2.7) 

The evolution for pj* was considered by Mueller in his original work . It obeys the linear BFKL 
equation POj- The n-dipole densities with n > 1, are necessary generated through the evolution 
even if the projectile is a single dipole. For this specific choice of the projectile, our linear equation, 
say, for can be identically reformulated as a nonlinear with respect to pf. This form of the 
equation for pf was studied in Refs. PHI]. We would like to stress, however, that Eq. ()2.7|) is 
more general, as it does not depend on any specific choice of projectile. Our approach cannot be 
rigorously applied to a hadron or nucleus as a projectile, as the latter need to account for high 
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dipole densities, in spite of this we view our finding as a first important step towards a symmetric 
description of target and projectile. 

So far we have discussed the evolution of the projectile only. Our goal is to compute the total 
scattering amplitude. Following Ref. ,3j the amplitude N is defined 

oo „ n 

NiY) = - J 7n(ri,6i...,r„,6„;Fo)p^(ri,6i...,r„,6„;r-ro) l[d\d\. 

n=l i=l 

(2.8) 

The amphtude for simultaneous scattering of n dipoles off the target is denoted by 7„. It has to 
be specified at the lowest rapidity (Yq). 7^ can be expressed through p^, the dipole densities in 
the target: 

7„(ri, 61,... ,r„,6„;yo) = 

(2.9) 

00 „ 

/ a'""(ri, 61, . . .r„, 6„ | fi, 61, . . . , f^, 6^) p^(fi, 61 . . . , 6„,; Fq) fid^hi . . . fmd^hm, 

m=l 

with 0"" '" being the properly normalised amplitude for n dipoles scattering off m dipoles, at the 
two gluon exchange level for each pair of dipoles. 

While the generating functional Z contains information about the projectile only, the target 
enters the amplitude through the functions 7^ (Pm)- Substituting ()2.9|) into ()2.8p we obtain: 

00 „ 

N{Y) = - i-^T J pl{fi,h...,fr„bn,Yo)ff^{n,h...,rr„b^;Y-Yo) (2.10) 

n, m=l 

n m 

X (j"™(ri, 61, . . . r„, 6„ I fi, 61, . . . , f^, &m) Y[ d'^Tid^hi W d^fjd'^hj. 

i=i j=i 

Eq. ()2.10|) appears symmetric with respect to the exchange of the target and projectile. In fact this 
symmetry is illusive. As soon as we require to obey the evolution given by Eq. ()2.7|) the boost 
invariance would impose a different evolution law (which we will not be able to write explicitly) 
for p^. We would like to emphasise the importance of the summation over n. It guarantees 
the unitarity, and boost invariance of the amplitude. It is challenging to include simultaneously 
high dipole density effects, in both the evolution of target and projectile. Eq. ()2.10|) would then 
account for all the effects under discussion in Ref. [13] . 

If ^0 = Eq. ()2.8|) yields the expression for the amplitude in which the whole evolution is 
through the projectile wavefunction. The dipoles produced via the evolution then multiple rescat- 
ter off the target resulting in the unitarity respecting amplitude. The unitarisation is achieved 
without saturation, as no high density effects are accounted by the evolution of the projectile. 

Alternatively we can take Yq = Y . This would bring the whole evolution into the target 
averaged amplitudes. By requirying the total amplitude to be independent of our choice of 
frame we will be able to deduce the evolution law for 7n(F)- As can be expected 7^ happens to 
obey the Balitsky's chain of hierarchy equations j3] projected onto dipole operators. 
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We now derive the evolution of 7„. We note that 
We obtain the hnear functional equation for the amplitude A^: 



dN CO n f 6 N 

— ^T7 = J2 Yl £rid^k 'yniri,bi,. . . ,rn,bn;Yo) [- ^ uj{ri) - — 

OisCI J- n=l i=l \ i=l " 



The equation (j2.1H) can be solved using the ansatz: N{Y, [7]) = A^(7i(F), 72(1^) . . .) such that 



dY S',„ dV i\ 



As all 7„ are independent of each other, Eq. ()2.11|) is equivalent to Balitsky's hierarchy (restricted 
to dipole operators): 

= - 2^ (^{ri) 7n(n, Oi,--- ,rn, bn) + 



E 

i=l 



d? r' 



2 TT {vi — r') 



/^2 



2 7„(ri, 61,... r', k + . . . ,r„, 6„) (2.12) 



- 7n+i(n, &1, . . . r , 6i H ^ , . . . [Ti - r ), b, + - 

The initial conditions for the evolution chain ()2.12j) are given by a target model which is supposed 
to provide 7„(F = Yq). All 7„ are defined as scattering amplitudes and hence are target averaged 
quantities. In principal, Eq. ()2.12|) governs the evolution of the target densities p^. Though we 
are not able to write an explicit equation for p^, it apparently differs from the evolution of p^ (Eq. 

The amplitude A^ can be rewritten in the following form 

00 „ n 

N{Y) = - E(-iry 7n(ri, 6i,...,r„, 6„;r-ro)p^(ri, &1, ...r„, 6„;yo) U n dH, . 

n=l i=l 

The Y evolution has been transferred from the evolution in the projectile wavefunction (p*^) to 
the evolution of the target averaged amplitudes (7^). Note that by construction N(Y) is boost 
or frame invariant. In other words, Eq. (j2.12j) can be obtained just by requirying N{Y) to be 
Yq independent. The functions p^{Y = Yq) have to be provided by the initial dipole distribution 
in the projectile. So far, our treatment was general. Note, however, the asymmetry between the 
target and projectile, as for the latter no high density effects are included. This asymmetry is 
reflected by the fact that p^ and 7„ obey different hierarchy equations. 
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We would like to emphasise that so far our discussion was quite general and valid for any 
projectile. Our results can be trusted until we violate our central assumption, namely the absence 
of nonlinear effects in the projectile wave function. 

We now make a choice of a projectile. Consider the canonical case of a dipole of the size r 
scattering off a generic target at the impact parameter b. Then the initial condition in the projectile 
wave function at Y = Yq is = P„=i = 5^(r — ri) 5^(6 — bi) while p^yi = Pn>i = 0. This 
choice of the projectile initial conditions implies 

Z{Y = Yo, r, b; [u]) = u{r,b) . (2.13) 

For the total amplitude we immediately obtain 

N{Y, r, b) = 71 (r, b; Y) 

For a generic target we cannot write down a closed form equation for A^. It was shown in Ref. 
[2] that for this particular choice of the projectile, Eq. ()2.4|) can be rewritten in the nonfunctional 
but non-linear form, reproducing the same equation for Z as in Ref. [T]: 

agOY 



Z\Y-r\b+ N Z\ Y- (r - r'), b - -; [u] 



The linear functional equation ()2.4j) is more general, however, as it is valid for any projectile. 

3 Target correlations 

Target correlations are defined by having 

7„(ri, 61, . . . r„, bn] Y = Yo) = C„(ri, 61 . . . r„, 6„) 7(ri, 61) . . . 7(r„, 6„) 

The coefficients C„ are ra-dipole correlation parameters. They can be also viewed as effective 
measures of target fluctuations. If we assume C„ be pure numbers independent of coordinates, 
then Eq. 1)2.111) is reduced to the following functional equation j2] 

^^^^Py = - J d'r' V,^,ir', [7(r')]) N{Y; [7]) (3.15) 

+ A / d'r'd'r" 7(r") , f N{Y; [7]) 

2 IT J ^'/(^j." _ j^'y d'-f{rj 

+ f(^J d'f^if) J d'r" d'r' V^^2{r', r", [j{r')]) N{Y- [7]) 



7 



Here F{n) = Cn/Cn-i- If we further simplify the correlations assuming F is constant 

a = Ci (3.16) 

then Eq. ()3.15j) can be easily reduced to a non-functional non-linear equation (ignoring the b 
dependence) : 

dN(r,Y) f <Pr' 

- "5 X / -- -— — X (3.17) 



dY JpP 2tt r'2 (r - r')^ 

[2N{r', Y) - N{r, Y) - F/Ci N{r', Y) N{r - r', Y)] , 

The initial condition is N(Y = Yq) = Ci 7. Note that unless target properties are specified, F 
and Ci are arbitrary numbers {F > Ci). 

Our solution first obtained in Ref. |2|, proves that under condition ()3.16|) the hierarchy ()2.12j) 
shrinks to a single nonlinear equation ()3.17|) . An identical result has been obtained recently in 
Ref. P by a direct analysis of the hierarchy equations, whereas the result F/Ci = 2 of Ref. [12] 
is a particular solution. 

The new aspect which we find as most important, is the relation between the target correlations 
and Balitsky's hierarchy, something which we were not aware of when preparing Ref. j2]. Eq. 
(|3.17p which we derived was viewed as a phenomenologically motivated ad hoc modification of the 
BK equation. By establishing in this letter a relation between the correlations and hierarchy, we 
attempt to correct the wrong attitude regarding our original work. 

If we consider dipole interaction as fully uncorrelated F = Ci = 1 then Eq. (j3.17p reduces to 
the BK equation. Its solution can also be written in the following form 

00 n 

iV(y, r, 6) = - 5](-l)"p^(ri, fei, ...r„, 6„;F-ro) 11^(^0,^,6.) ^r^dH,. 

71=1 i=l 

So far the approximation F = const looks pure mathematical. In Ref. [2], a realistic model 
for target correlations was proposed. In that model the target was considered to be a nucleus (not 
necessary very heavy) but the approach can be viewed more generally. The target correlations 
were estimated based on pure counting arguments independent of the dipole coordinates. The 
model of Ref [2] leads to a constant F, providing a physically intuitive realisation of a more 
formal structure outlined above. 
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